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We analyze the macroscopic dynamics of a Bose gas axially confined in an optical lattice with a 
superimposed harmonic trap, taking into account weak tunneling effect. Our results show that upon 
transition to the quasi-two-dimensional (2D) regime of the trapped gas, the 3D equation of state 
and equilibrium density profile acquire corrections from 2D many-body effects. The corresponding 
frequency shift in the transverse breathing mode is accessible within current facilities, suggesting a 
direct observation of dimensional effects. Comparisons with other relevant effects are also presented. 
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The interplay between dimensionality and quantum 
fluctuations in two-dimensional (2D) strongly correlated 
quantum systems [1] has long been recognized to give rise 
to remarkable phenomena like high- Tc superconductivity 
[2] and the long-sought Berezinskii-Kosterlitz-Thouless 
(BKT) transition [sj. Recent extraordinary realization 
of quasi-2D ultracold Bose gases [4], where tightly con- 
fined axial kinematics manifests as 2D features in pair 
collisions fsj, has cast new light in understanding low- 
dimensional many-body systems [6[. 

In view of the important role played by dimension- 
ality, their spectroscopic diagnostics is highly interest- 
ing [7|. The interest particularly originates from the 
sensitivity of collective frequencies to equation of state 
which henceforth establishes their measurements as pre- 
cise tests to the many-body physics. For example, the 
universal breathing mode of a 2D harmonically trapped 
Bose gas with ^^^(r) interaction, first noticed via the 
Castin-Dum-Kagan-Surkov-Shlyapnikov (CDKSS) scal- 
ing ansatz [8|, |9[, reveals the hidden Pitaevskii-Rosch 
symmetry (PRS) in the associated classical field 

theory. Whereas, the frequency shift in this mode 
can provide a signature of quantum anomaly emerging 
upon quantization when quantum fluctuations signifi- 
cantly modify the scattering length [12]. 

In this Letter, we are inspired to discuss visualization 
of the effect of dimensionality in collective excitations 
of an experimentally favored quasi-2D optically trapped 
Bose gas (seen in Fig. 1), taking the quantum tunnel- 
ing into account. Collective oscillations in the presence 

of optical lattice has been previously studied intensively 

1-2 2-1 

for the effective Hamiltonian = + 2^ + 

f Ej H + ^1) + ^^'] +^ E,<fe where ~g is 

the lattice-renormalized 3D coupling constant and m* 
is the effective mass [l3|, Q. The corresponding su- 
perfluid hydrodynamic analysis predicted that ax- 
ial optical potential mark no effect on transverse modes 
for an elongated trap, which was later experimentally 
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FIG. 1. (color online). Schematic picture of an optically 
trapped quasi- 2D Bose gas. A ID optical lattice is along the 
horizontal axis (z axis) with a BEC in an elongated harmonic 
trap, with axial (radial) frequency Uz (<^±). The BEC is thus 
confined to an array of narrow potential pancakes. 



tested by Fort et al [15]. However, this scenario will ex- 
perience fundamental modifications when transiting to 
the quasi-2D regime. As we shall show in this Letter, 
upon the onset of "frozen" axial kinematics in tight op- 
tical lattice, a 2D-peculiar position-dependent correction 
Hi = ^j^j. gi (rj) 5^ i'^jk) emerges from the perspective 
of many-body theories [5]. Its impression on collective 
frequencies can be foreseen, for example, from the equa- 
tion of motion 
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for the excitation operator Fb = +^2^) trans- 

verse breathing mode. Equation ([T]), being equally valid 
classically as it is quantum mechanically, immediately 
leads to following statements: (i) for a pure 2D system 
within the classical field description, only Hq survives on 
the right side of Eq. ([T]) and the mode operator oscillates 
universally with 2uj±, as required by PRS [lo|; (ii) this 
universal oscillation with 2uj± also persists in a very elon- 
gated 3D dilute Bose gas described by Hq {uJz/cjO± <C 1), 
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as predicted by Ref. [13|; (iii) the oscillation frequency 
can be shifted from 2uj±_ by the emerging correction ^i(r). 

In what follows, we analytically calculate the col- 
lective excitations of a quasi-2D Bose gas tightly con- 
fined by an optical lattice Vopt = sEiis\v?{qBz)^ with 
a superimposed cylindrically symmetric harmonic trap 
Vho (r) = ^{uj\_x^ + ^1^^ + ^1^'^)-> as shown in Fig. 1. 
The lattice period is fixed hy qb = tt / d with d being the 
lattice spacing, s is a dimensionless factor labeled by the 
intensity of a laser beam and Er = h^q^/2m is the recoil 
energy with Hqe being the Bragg momentum. 

Our starting point is the linearized hydrodynamic 
equation for density fiuctuations Sn{r^t)^ generalized 
straightforwardly from Ref. [13] to the quasi-2D regime. 
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with V = ^Vx, V^^/m/m*^ . The key ingredient of Eq. 

([2]) is the zero-temperature local chemical potential ij^q2d 
of the quasi-2D Bose gas under consideration. Equation 
([2]) is justified by sufficiently weak tunneling which is 
nevertheless nonnegligible to ensure full coherence of the 
order parameter between different wells [13], and by as- 
suming the Thomas-Fermi (TF) limit and local density 
approximation [7]. Here, the 3D density n (r) is deter- 
mined from /io = l-i^Q2D [^(r)] + Vho (r), where /io is the 
ground state value of the chemical potential, fixed by the 
proper normalization of n (r) . 

At the core of hydrodynamic analysis on collective os- 
cillations is the knowledge of the equation of state. In 
order to determine jiiQ2D in Eq. ([2]), we start from the 
grand partition function ^ Z = J D [t/j'^.tp] g-'^t^*'^]/^ 
of an optically trapped quasi-2D Bose gas in the absence 
of harmonic potential, where 



S[il)\^l)]= j dr j d^Y^I)''{Y,T) 



+ Kpt(r) + ||V^(r,r) 



dr 2m 
i^{v,T) (3) 



is the action functional of [?/^* (r, r) , (r, r)] which collec- 
tively denote the complex functions of space and imagi- 
nary time r. Here, Qe abstractly stands for the two-body 
coupling constant in an axial optical confinement. Us- 
ing the path-integral approach [16|, one finds within the 
tight-binding approximation and Bogoliubov theoretical 
framework [l7| the ground state energy Eg^ 
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Here, t denotes the tunneling rate. 



n refers to the 
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condensate density, and ge = Qe 

g(.(dl \/2j\g\ where w(u) = exp 
is a variational Gaussian anstaz and the ratio d/a 
7r5^/^ exp (—l/4:\/s) minimizes the free energy functional 
with respect to a |l4|. 

The ground state energy in Eq. (j4]) experiences a 
lattice-induced dimensional crossover governed by the pa- 
rameter 2t/gen. In the limit 2t/gen ^ 1, one finds the 
system exhibiting anisotropic 3D behavior and F{x) 
32/15 v^. Whereas, the extreme 2t/gen <C 1 corresponds 
to the quasi-2D regime where axial atomic motion is 
frozen to zero-point oscillations and = ^ — ^ logx 

is approached exactly. We stress that above analysis is 
justified by the weak but nonegligible tunneling effect 
guaranteed by 1/A^^ < t/geU with Nt being the number 
of atoms per optical well. 

Upon transiting from 3D to quasi-2D regime 1/Nt < 
2t/gen <C 1, the axial optical lattice imprints its effect 
in binary atomic collisions via strongly restricting axial 
kinematics that manifests as a 2D character of the rela- 
tive motion of colliding atoms at large separation. This 
gives rise to an effective coupling constant sensitive to the 
lattice parameter which is approximately given by [HI, [HI 

^ 2^h^d 1 

m a2D/a3D + (1/V27r) In [l/n2Da^^] ' 

where n2D = ^d is the surface density and a2D = cr is the 
effective 2D scattering length. The logarithmic density- 
dependent term in Eq. (|6j) is typical of 2D many-body 
effects, whose relative importance is governed by the ra- 
tio a2D/ci?,D which henceforth controls the dimensional 
crossover in hard-core interactions [5]. For a2Dl olzd ^ 1, 
one finds pure 3D collision and Eq. ([6j) converges to the 
lattice-renormalized density independent coupling con- 
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Whereas in the opposite limit 



ci2D/oid,D <C 1, 2D scattering dominates which leads to 
9e = g2Dd with g2D = (47r/i^/m)/ln (l/n2Dtt2D) being 
the density-dependent effective coupling constant pecu- 
liar to pure 2D Bose gases. Hereafter, we focus on what 
we shall call the weak quasi-2D regime where the small 
parameter a2D/ci^D measures the deviation from pure 3D 
scattering due to emerging 2D features. Linearly expand- 
ing Eq. ([6]) with respect to a^D/ci2D^ one obtains 
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where the term showing the logarithmic dependence on 
the gas parameter constitutes the leading 2D correction 
to 3D coupling constant g. 

The jiiQ2D can now be readily determined via /i = 
dEg/dN from Eq. (|4]), together with proper asymp- 
totic analysis. We first note that in the asymptotic 
3D regime where 2t/gn ^ 1 and a^D/ci2D ^ 1, our 
analytical solution (i = gn 1 + (32m*/3v^m) Y^ag^n 

is consistent with the 3D Lee- Huang- Yang (LHY) re- 
sult [19]; whereas, in the opposite pure 2D limit where 
2t/gn <C 1 and a^D/ci2D ^ 1, our asymptotical re- 
sult for the chemical potential of a 2D Bose gas fi = 

4.7vh'^n2D/m L _ ln(ln(l/n2i?a2i^))-B ^ ln(ln(l/n2ija2i^))-B 
\lnn2Dal^\ [ In(l/n2i?a2^) In'^ {l/n2Dal^) 

with B = l — ln {mt / n2D27ih^) stands in good agreement 
with Ref. [ioj. Thence, after applying similar schemes to 
the weak quasi-2D regime where 2t/gn<^ 1 and a^D/ci2D 
is small, one finds 
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Rewriting /iQ2D = [1 + (^)], we thus identify 
k2D(n) = asD/{V27ra2D)[l/2-\n {l/n2Dal^)] as the 
first correction to the 3D mean-field (MF) equation of 
state arising from the 2D effect. 

From Eq. ([8|), the equation for the 3D ground state 
density can be solved by iteration yielding 
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with Utf (r) = (/io — Vext (r)) /g being the 3D TF den- 
sity. Eq. clearly shows that, because of the weak 
coupling between adjacent wells, the 2D corrections in 
the local chemical potential in Eq. ([8]) are transferred to 
the 3D stationary shape of cloud. 

Substituting Eqs. and into Eq. ^ and only 
retaining terms linear in /c2d(^), we obtain 
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Equation (fTQ|) in the absence of k2D is just the familiar 
3D hydro dynamic equation in the presence of ID opti- 
cal lattice [13|. Against this background, the addition of 
terms on the right side of Eq. (p!Q|) presents a pertur- 
bation. The resulting fractional frequency shift, to the 
leading order, is given by 
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where integrals extend to the region where ut f is positive 
[2T] . An important feature in Eq. (pT]) is the dependence 
of 5uj/uj on the derivative dk2D/dn rather than /c2d(^)- 



The consequence is that the leading order correction aris- 
ing from 2D effect to 3D MF collective frequency shows 
no logarithmic density dependence. 

According to Eq. (pT]) . the surface modes that sat- 
isfy V^^n = are not perturbed by the 2D effect in 
the vicinity of 3D regime. Hence, in order to observe 
dimensional effects, one has to focus on small com- 
pressional oscillations. Our primary mode of interest 
is the transverse breathing mode in a very elongated 
trap {^/mJmF uo z / ^ 1. <C 1). Substitutions of 5n{v) ^ 
— with Rtf = 2pL^ j muj\ being the trans- 

verse TF radius and uo = 2uj^ into Eq. (fTTj) yield the 
fractional shift 

buo 
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Equations (pTj) and (fT2]) consist of the major results 
of this paper. In typical experiments to date [22[, the 
relevant parameters are given by 3D scattering length 
<^3D = 5.31nm and the lattice period d = 297. 3nm. The 
frequency shift in Eq. ([12]) can be reached ~ 0.48% for 
s = 4. Given an accuracy of ~ 0.3 — 0.4% in measur- 
ing collective frequencies within current facilities [23|, 
the 2D correction to the transverse breathing mode is 
well in reach in relevant experiment conditions. More- 
over, this effect can be enhanced via adjusting lattice 
parameter and using Feshbach resonance. We have also 
taken a look at the lowest compression mode in a disk- 
like geometry {^/m/m*uOz/uo± ^1), which is along the 
axial direction with the zeroth order dispersion given by 
uo = ^y/vfij m* ) \f?iUJz and density oscillations of the form 
(5n(r) — Z\plZ^ where Ztf = \/ 2m* ji^ / muol is the 
TF radius along the axial direction. Straightforward cal- 
culations yield — = 1 223_^3P showing an amplified 

^ w 6V27r rn a2D ' ^ ^ 

2D effect due to the increased inertia along the direction 
of the laser. 

The frequency shift in Eq. ([T2]) should be compared 
with other relevant corrections in actual experiments, like 
finite size, nonlinearity, thermal effects and vortex. Fi- 
nite size effects originate from kinetic energy pressure 
typically ignored in the TF scheme. Its consequence on 
the transverse breathing mode can be analyzed from Eq. 
(dJ), averaged over the many-body wavefunctions in ac- 
cordance with the experimental measurements. For very 
elongated geometry where nearly spatial invariance in ax- 
ial direction implies Ekinz = 0, calculations by GP theory 
directly give uo = 2uo±^ showing no finite size effect. In 
general case where uOz is finite and Ekinz 7^ 0, a sum rule 
approach within single-mode approximation gives uo = 
Y^ms/mi with rns = (Sfi^/m^) {Ekin± Eho± ^ Eint) 
and mi = {2h^ /m)N {Fb) respectively being the cubic 
energy weighted and the energy weighted moments of the 
dynamic structure factor nip = | (0|Fb |n) p(flCJno)^• 
Using the viral identity Ekin± — Eho± + Eint = for the 
ground state in GP description, one again finds uj = 2uj± 
unaffected by the finite size effect. 
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In actual experiments where the amphtude of os- 
cihation can not be arbitrarily small, effects of non- 
linearity can shift the collective frequency approx- 
imately hy 5uj/uj A^5 Here A is the 
fractional oscillating amplitude of the harmonically 
confined atomic cloud which can be tuned to less 
than 10%. In our case, the coefficient 5 is cal- 
culated as 5 = t ^?7-^^^.^7'^%-,^^ 



is supported by the NSF of China (Grant No. 11004200) 
and IMR SYNL-TS Ke Research Grant. 



15 



16 (g_ 
3 g+ 



-3 r 



20 



q+-q- 



-10AV + 37A2 



9A2- 



-| I 15 A 



(g--4) 



4 

11(7+ - 54] with q± = 2 



|A^ T ^V9A^ - 16A2 + 16, which asymptotically var- 
nishes in the limit A = uJz/^± 0. The effect of nonlin- 
earity is therefore very small. 

The consequence of thermal effects on collective ex- 
citations has been experimentally investigated in Ref. 
[11]. The observed unusually small damping rates and 
frequency shifts in the transverse breathing mode of an 
elongated condensate have been attributed to an acci- 
dental degeneracy between the condensate and thermal 
cloud oscillation frequencies [25]. This accidental sup- 
pression of Landau damping may offer an experimental 
control of the thermal effect on the frequency shift. 

The issue of vortex is closely related to operational 
anisotropy in excitation schemes. In actual exper- 
iments, the breathing mode is generally excited via 
quenching the transverse harmonic trap frequencies in 
phase by SuJx(^y^ = 5uo <C uo^. However, this scheme 
will be perturbed by an out-of-phase operation with 
5uJx = and 5ujy = duo , and the resulting excita- 
tion operator F = muj±{Suj -\- 5uj )/2 x Fb + muj±{5uj — 
5uj )/2 rf (F2,2 + ^2,-2) iXim being spherical harmon- 
ics) gives rise to additional excitations of quadrupole 
modes. Such situation is further spoiled by the presence 
of quantum vortex which splits the m = ±2 quadrupole 

modes by approximately cj+-a;_ = 

We note that this split becomes infinite when A ^ 0, sug- 
gesting an important point of observing our result in Eq. 
(p!2]) is to avoid the anisotropy in perturbation schemes 
and presence of a vortex. 

In conclusion, our results show that 2D many-body ef- 
fects can be visible in the frequency shift of the transverse 
breathing mode for an optically trapped Bose gas transit- 
ing from 3D to quasi-2D regime. Observing dimensional 
effect directly would present an important achievement in 
revealing the interplay between dimensionality and quan- 
tum fluctuations in low-dimensional strongly correlated 
quantum systems associated with BEG. 

Noting that compared to the fruitful work on a BEG 
along the 3D, quasi- ID and ID dimensional crossover 
[l^l, much less have been reported in the transition re- 
gion from the 3D to quasi-2D and 2D regimes where many 
unique phenomena are known to arise. More comprehen- 
sive study along this line is highly desirable. 

We thank Biao Wu for helpful discussions. This work 
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